Approximate and exact nodes of fermionic wavefunctions: coordinate transformations 

and topologies 
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A study of fermion nodes for spin-polarized states of a few-electron ions and molecules with 
s,p,d one-particle orbitals is presented. We find exact nodes for some cases of two electron atomic 
and molecular states and also the first exact node for the three-electron atomic system in 4 S(p 3 ) 
state using appropriate coordinate maps and wavefunction symmetries. We analyze the cases of 
nodes for larger number of electrons in the Hartree-Fock approximation and for some cases we 
find transformations for projecting the high-dimensional node manifolds into 3D space. The node 
topologies and other properties are studied using these projections. We also propose a general 
coordinate transformation as an extension of Feynman-Cohen backflow coordinates to both simplify 
the nodal description and as a new variational freedom for quantum Monte Carlo trial wavefunctions. 
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I. 



INTRODUCTION 



The problem of fermion nodes is one of the most in- 
triguing challenges in quantum simulations of fermionic 
systems by stochastic methods such as quantum Monte 
Carlo (QMC)i£. In QMC the many-fermion wavefunc- 
tion is represented by an ensemble of sampling points 
(walkers) in the space of fermion coordinates. The walk- 
ers are propagated according to the matrix elements of 
the projector exp(— tH) where r is a real parameter and 
H is a Hamiltonian. It is straightforward to show that 
the propagated wavefunction solves the time-dependent 
Schrodinger equation in imaginary time r = it and con- 
verges to the ground state for r — > oo. Unfortunately, for 
fermions such a straightforward sampling process runs 
into difficulties and the projection becomes very inef- 
ficient. The ensembles of walkers which initially sam- 
ple negative and positive parts of the wavefunction are 
independent and asymptotically converge to the same 
(bosonic) distribution with an exponential growth of er- 
ror bars for fermionic observables. In electronic structure 
QMC calculations 3 this well-known fermion sign problem 
is circumvented by the fixed-node approximation which 
restricts the negative and positive walkers into separate 
regions of space defined by an approximate fermion node 
(zero boundary) of the best available trial/ variational 
wavefunction. That guarantees stability of the statistical 
error bars at the price of a fixed-node bias. The fixed- 
node bias is proportional to the square of the nodal dis- 
placement error and therefore in typical electronic struc- 
ture calculations the resulting bias is rather small. Even 
for hundreds of electrons, Hartree-Fock or multi-reference 
Hartree-Fock (HF) nodes lead to impressive accuracy 
when used within QMC. The fixed- node QMC calcula- 
tions typically provide about 95% of the correlation en- 
ergy in real systems^ such as molecules, clusters, solids, 
etc. 

At a more fundamental level, knowledge of the exact 
node enables one to eliminate the fixed-node bias com- 



pletely, and the exact energy can be calculated in time 
which scales as a low-order polynomial in the number of 
particles. Therefore elimination of the fixed-node error 
remains one of the intriguing possibilities for employing 
QMC to attack a number of important many-body prob- 
lems which require accuracy beyond a few-determinant 
Hartree-Fock nodes. 

Let us assume a system of spin-polarized electrons 
described by a real wavefunction ^f(R) where R de- 
notes the electron spatial coordinates. The exchange of 
an electron pair with labels denoted as Pij, gives 
^(i?) = —ty (P^ R) . Consequently, the antisymmetry im- 
plies that there exists a subset of electron configurations, 
called a fermion node, for which the wavefunction is zero. 
Let us eliminate the regions in which the wavefunction 
vanishes because of other reasons (eg, external potential); 
then the fermion node is given by an implicit equation 
>F(i?) = 0. In general, the fermion node is a (ND-1)- 
dimensional manifold (hypersurface) assuming that we 
have N fermions in a D— dimensional space. The fermion 
nodes of small systems, mostly atoms, were investigated 
in several previously published papers^SiLS. The gen- 
eral properties of fermion nodes were analyzed in an ex- 
tensive study by Ceperley 10 which included a proof of the 
tiling property and generalizations of the fermion nodes 
to density matrices. We mention two of the results which 
will be used later. 

i) Tiling property for the nondegenerate ground state: 
Let us define a nodal cell £l(i?o) as a subset of configura- 
tions which can be reached from the point Rq by a contin- 
uous path without crossing the node. The tiling property 
says that by applying all possible particle permutations 
to an arbitrary nodal cell of a ground state wavefunction 
one covers the complete configuration space. Note that 
this does not specify how many nodal cells are there. Fur- 
thermore, symmetry of the state is also symmetry of the 
node and tiling property is valid for any non-degenerate 
ground state within the given discrete symmetry. 

ii) If two nodal surfaces cross each other they are or- 
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thogonal at the crossing. If n nodal surfaces cross each 
other, the crossing angles are all equal to n/n. 

In addition, it was numerically shownltt that the nodal 
cells for a number of finite free particle systems are max- 
imal, ie, all regions with the same sign of wavefunc- 
tion are interconnected. The fermion nodes for degen- 
erate and excited states were further studied by Foulkes 
and co-workersii. Recent interesting work by Bres- 
sanini, Reynolds and Ceperley revealed differences in the 
nodal surface topology between Hartree-Fock and corre- 
lated wavefunctions for the Be atom explaining the large 
impact of the 2s, 2p near-degeneracy on the fixed-node 
QMC energy-^. 

The remainder of this paper is organized as follows: in 
Sec. [n]we discover new exact fermion nodes for two and 
three-electron spin polarized systems. In Sec. IHII we cat- 
egorize the nodal surfaces for the several half-filled sub- 
shells relevant for atomic and molecular states. In Sec. 
II VI we suggest a general particle position transformation, 
both as a tool to simplify the description of the nodes 
and also as a possible new variational freedom for trial 
wavefunctions. Finally, in the last section we present our 
conclusions and suggestions for future work. 



II. EXACT NODAL SURFACES 

We assume the usual electron-ion Hamiltonian and we 
first investigate a few-electron ions focusing on fermion 
nodes for subshclls of one-particle states with s,p,d, f... 
symmetries using variable transformations, symmetry 
operations and explicit expressions for the nodes. 



A. Three-electron quartet 4 5'(p 3 ) state 

Let us first analyze a special case with r% — r 2 and 
r 23 — r 3i- R is then easy to see that the inversion around 
the origin with subsequent rotations is equivalent to the 
exchange of two particles, say, 1 and 2 (FigQ}. Therefore 
for this particular configurations of particles the combi- 
nation of parity and rotations is closely related to the 
exchange symmetry. The illustration also shows that the 
six distances do not specify the relative positions of the 
three electrons unambiguously. For a given set of the dis- 
tances there are two distinct positions, say, of the elec- 
tron 3, relative to the fixed positions of electrons 1 and 2 
(see Fig0 and compare positions 3 and 3" of the third 
electron. 

In order to analyze the wavefunction in an unambigu- 
ous manner it is convenient to define new coordinates. 
Let us denote r^" 2 = t\ + r2,?i2 = |r]^|, together with 
the customary r 12 = i*i — r 2 ,r 12 = |r 12 | . We can now 
introduce the following map of the Cartesian coordinates 

(ri,r 2 ,r 3 ) -> (rf 2 , r 12 , r 3 , cos a, cos/3, 7, fi) (1) 




FIG. 1: Inversion and two subsequent rotations of three par- 
ticles: (a) Original and inverted (primed) positions; (b) Posi- 
tions after the rotation Ri in the plane given by the particles 
1,2 and the origin; (c) Positions after the second rotation R2 
around the ri + r2 axis. Note that the original positions of 
the particles 1 and 2 are exchanged. 

with definitions: cos a = r 3 • (ri x r 2 )/(r3|ri x r 2 |), 
cos/3 = • ri 2 /(r 1 " 2 ri 2 ) and 7 being an azimuthal an- 
gle of r 3 in the relative coordinate system with unit 
vectors = r+/r+, e z = n x r 2 /\r 1 x r 2 |, e y = 
e z x e x . For completeness, f2 denotes three Euler an- 
gles which fix the orientation the three-particle system 
in the original coordinates (eg, two spherical angles of 
ri x r 2 and an azimuthal angle of rj" 2 ). Since the an- 
gles fi are irrelevant in S symmetry, the first six vari- 
ables fully specify the relative positions of the three 
particles and the wavefunction dependence simplifies to 
^( r i2> r i2i r 3i cos a : cos/3, 7). Consider now two symme- 
try operations which change the sign of the wavefunc- 
tion and keep the distances unchanged: parity Pj and 
exchange Pi 2 between particles 1 and 2. The exchange 
flips the sign of all three cos a, cos /3, 7 while the parity 
changes only the sign of cos a. The action of P/Pi 2 on W 
leads to 

*&(..., cos a, — cos/3, —7) = \E , (..., cos a, cos/3, 7) (2) 

showing that the wavefunction is even in the simultane- 
ous sign flip (cos/3,7) — > (— cos/3, — 7). Applying the 
exchange operator Pi 2 to the wavefunction and taking 
advantage of the previous property gives us 

— cos a, cos /3, 7) = —"&(..., cos a, cos/3, 7) (3) 

suggesting that there is a node determined by the condi- 
tion cos a = 0. It is also clear that the same arguments 
can be repeated with exchanged particle labels 2^3 and 
3 <-> 1 and we end up with the the same nodal condition: 
r 3 • (iq x r 2 ) = 0. This shows that the node is encountered 
when all three electrons lie on a plane passing through the 
origin. Now we need to prove that this is the only node 
since there might possibly be other nodal surfaces not 
revealed by the parametrization above. The node given 
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above clearly fulfills the tiling property and all symme- 
tries of the state. Furthermore, the state is the lowest 
quartet of S symmetry and odd parity (lower quartets 
such as Is2s3s,ls2s2p, and \s2p 2 have either different 
parity or symmetry) and for the ground state we expect 
that the number of nodal cells will be minimal. This 
is indeed true since the node above specifies only two 
nodal cells (one positive, one negative): an electron is 
either on one or the other side of the nodal plane passing 
through the remaining two electrons. Futhcrmore, any 
distortion of the node from the plane necessarily leads to 
additional nodal cells (see Fig. |2J which can only increase 
energy by imposing higher curvature (kinetic energy) on 
the wavefunction. This is basically the Feynman's argu- 




FIG. 2: (a) An illustration of an artificial distortion of the 
planar ground state node for the quartet state (see text); (b) 
The original and parity transformed distorted node; (c) Fi- 
nally, a subsequent rotation of the inverted distortion nec- 
essarily leads to a new nodal pocket which is artificial for 
the ground state. In fact, nodes with similar topologies are 
present in excited states—. 

ment from the proof demonstrating that the energy of 
fermionic ground state is always above the energy of the 
bosonic ground state (and also essentially the same ar- 
gument as used for the proof of the tiling property"). In 
fact, as we show in another paper, higher excited states 
of this symmetry have additional nodes, as expected^. 
Given all the arguments above we conclude the proof 
that the plane is the exact node. Note that it is identical 
to the node of Hartree-Fock wavefunction of 2p orbitals 
given by ^hf — det{p(r)x, p{r)y, p(r)z} where p(r) is a 
radial function. 

The coordinate transformation above is not the only 
one which can be used to analyze this state. The high 
symmetry of the problem enables us to find an alterna- 
tive coordinate map with definitions of cos (3 modified to 
cos^' = [(n x r 2 ) x rf 2 ] • r 12 /(|(ri x r 2 ) x r^ 2 ||r^ 2 |) and 
7 to 7' by redefinition of e z to e' z = [(iq x r 2 ) x iq 2 ] x 
r i2]/l( r i x r 2) x r i2] x r^JI and e' y = e' 2 ■ e x . In the 
redefined coordinates the search for the node simplifies 
to an action of Pn on ^(rjtj, ri 2 , r$, cos a, cos (3' , 7') 

*(..., -cos a, cos 0', 7') = -*(..., cos a, cos /5 ', 7') (4) 

since the distances and cos /?', 7' are invariant to P\i . Ob- 
viously, this leads to the same nodal condition as derived 
above. 



It is quite interesting to compare these two coordinate 
maps with /3, 7 and 7'. Although parity and exchange 
are independent operators the analysis above shows that 
in an appropriate coordinate system they imply the same 
nodal surface. Both these operators cause an identical 
sign change of the wavefunction indicating thus a spe- 
cial symmetry of the i S(p 3 ) ground state node which is 
higher than would be expected solely from antisymmetry. 
Similar observation was made in a study of fermion node 
in another case of two electron atomic stat oTi 14 . 

B. Two-electron triplet A P{p 2 ) and 3 E s (7r 2 ) states 

Apparently, the exact node of this case was derived in 
a different context by Breit in 193Cfei^. Here we offer an 
independent proof which enables us to apply the analy- 
sis to some molecular states with the same symmetries. 
The exact node for the 3 P(p 2 ) state can be found in a 
similar way as in the case of quartet above. The state 
has even parity, cylindric symmetry, say, around z-axis, 
and is odd under rotation by ir around x,y axes, R(irx), 
R(-ny). The mapping of Cartesian coordinates which en- 
ables to analyze the wavefunction symmetries is given by 

(ri,r 2 ) (r+ 2 , r 12 , cos uj, cos (3, ip, tp r ) (5) 

where cosw = zo ■ (rq x r 2 )/|iq x r 2 | with Zo being the 
unit vector in the z-direction and tp' being the azimuthal 
angle of iq x r 2 ; <p' can be omitted due to the cylindric 
symmetry. Further, tp is the azimuthal angle of r^ 2 hi 
the relative coordinate system with the x-axis unit vec- 
tor given by a projection of zo into the plane defined 
by iq,r 2 , ie, e x = z 0p /\z 0p \, e z = (iq x r 2 )/|iq x r 2 | 
and e y — e z x e x . Action of PiP^R^x) reveals that 
the wavefunction is invariant in the simultaneous change 
(cos f3,tp) — > (— cos (3, — ip). This property and action of 
P12 to the wavefunction together lead to 

*(..., -cosw,...) = -*(..., cosw, ...) (6) 

with the rest of the variables unchanged. The node 
is therefore given by cosw = and is encountered 
when an electron hits the plane which contains the z- 
axis and the other electron. As in the previous case 
the nodal plane fulfills the tilling property and mani- 
festly divides the space into two nodal cells so that we 
can conclude that this node is exact. The exact node 
again agrees with the node of Hartree-Fock wavefunction 
^ = a\et{p(r)x, p(r)y} . 

The fixed-node QMC energies for the 4 S(p 3 ) and 
3 P(p 2 ) cases derived above were calculated for a nitrogen 
cation with valence electrons in these states. The core 
electrons were eliminated by pseudopotential 21 . The trial 
wavefunction was of the commonly used form with sin- 
gle HF determinant times a Jastrow correlation factoA 
Note that the pseudopotential nonlocal s— channel docs 
not couple to either odd parity S state or even parity 
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TABLE I: Total energies (a.u.) of N+, N++ and N+++ ions 
with core electrons eliminated by pseudopotentials. The en- 
ergies are calculated by variational (VMC) and fixed-node 
diffusion (DMC) quantum Monte Carlo and Configuration 
Interaction (CI) methods. The HF energies are given as a 
reference for estimation of the correlation energies. 



State 


HF 


CI 


VMC 


DMC 


3 P(P 2 ) 


-5.58528 


-5.59491 


-5.59491(2) 


-5.59496(3) 


"Sip 3 ) 


-7.24716 


-7.27566 


-7.27577(1) 


-7.27583(2) 


^(sp 3 ) 


-8.98570 


-9.02027 


-9.01819(4) 


-9.01962(5) 



P{p 2 ) state so that that the nonlocal contribution to the 
energy vanishes exactly. 

In order to compare the fixed-node QMC calculations 
with an independent method we have carried out also 
Configuration Interaction (CI) calculations with ccpV6Z 
basis2£ (with up to three g basis functions) which gen- 
erates more than 100 virtual orbitals in total. In the 
CI method the wavefunction is expanded in excited de- 
terminants and we have included all single, double and 
triple excitations. Since the doubles and triples include 
two- and three-particle correlations exactly, the accuracy 
of the CI results is limited only by the size of the basis 
set. By comparison with other two- and three-electron 
CI calculations we estimate that the order of magnitude 
of the basis set CI bias is w 0.01 mH (miliHartree) for two 
electrons and w 0.1 mH and for three electrons (despite 
the large number of virtuals the CI expansion converges 
relatively slowly^S in the maximum angular momentum 
of the basis functions, in our case l ma x — 4). The pseu- 
dopotentials we used were identical in both QMC and CI 
calculations. 

The first two rows of Tab[I] show the total energies 
of variational and fixed-node DMC calculations with the 
trial wavefunctions with HF nodes together with results 
from the CI calculations. For 3 P(p 2 ) the energies agree 
within a few hundredths of mH with the CI energy being 
slightly higher but within two standard deviations from 
the fixed-node QMC result. For i S(p 3 ) the CI energy 
is clearly above the fixed-node DMC by about 0.17 mH 
as expected due to the limited basis set size. In order 
to illustrate the effect of the fixed-node approximation 
in the case when the HF node is not exact we have also 
included calculations for four electron state 5 S(sp 3 ) (for 
further discussion of this Hartree-Fock node see part III. 
below). For this case, we estimate that the CI energy 
is above the exact one by w 0.3 mH so that the fixed- 
node energy is significantly higher than both CI and exact 
energies. Using these results we estimate that the fixed- 
node error is w 1 mH, ie, close to 3% of the correlation 
energy. 

Since in the p 2 case we have assumed cylindric symme- 
try, the derived node equation is applicable to any such 
potential, eg, equidistant homonuclear dimer, trimer, etc, 
with one-particle orbitals tt x , n y which couple into the 



triplet state 3 S s (7r a ;7r y ). 

Note that the parametrization given above automat- 
ically provides also one of the very few known exact 
nodes in atoms so far—, ie, the lowest triplet state of 
He 3 S(ls2s). The spherical symmetry makes angles ui 
and f irrelevant and simplifies the two-electron wave- 
function dependence to distances r\,T2,r\2 or, alterna- 
tively, to ri2, rjtj, cos/3. Applying P\ 2 to wavefunction 
v I'( r 'i2,fi2 7 C0S /3) leads to 

-y(r 12 ,rf 2 , cos 0) = *(ri 2 ,r+,-cos/3) (7) 

so that the node is given by the condition cos/3 = 0, ie, 
r\ - r 2 = 0. 

In addition, the presented analysis sheds some light 
on the He 3 P(ls2p) state node which was investigated 
before 7 - as having higher symmetry than implied by the 
wavefunction symmetries. The symmetry operations re- 
veal that the wavefunction depends on | cos w | and that 
the node is related to the simultaneous flips such as 
(cos/3, ip) — > (— cos/3, — if) or angle shifts f — ► f + 7r. 
Since, however, two of the variables are involved, the 
node has a more complicated shape as the previous study 
illustrates 7 . In order to test the accuracy of the HF node 
we have carried out a fixed-node diffusion Monte Carlo 
calculation of the He 3 P(ls2p) state. The resulting total 
energy of -2.13320(4) a.u. is in an excellent agreement 
with the estimated exact value of -2. 13316 s which shows 
that the HF node is very close to the exact onei^. 



III. APPROXIMATE HARTREE-FOCK NODES 

It is quite instructive to investigate the nodes of half- 
filled subshells of one-particle states with higher angular 
momentum. 




FIG. 3: (Color online) The 3D projected Hartree-Fock node 
of 6 S(d 5 ) state which is an elliptic cone (left and right pic- 
tures). The middle picture illustrates a case when two pairs 
of two electrons lie on orthogonal planes which pass through 
the origin. This two-plane node is of lower dimension because 
of the additional condition on positions of the electrons. It 
appears as a crossover between the cones with different ori- 
entation (left and right pictures). The small spheres show 
the positions of the four electrons while the line denotes the 
z— axis. 
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A. Approximate Hartree-Fock node of 6 5'(d 5 ) state 

The HF determinant wavefunction for 6 S(d 5 ) is given 
*fff = ^ = iP{ri)det{2z 2 - x 2 - y 2 , x 2 - y 2 , xz, yz, xy} 

where p{ri) is the radial part of the d-orbital and we 
assume that all the orbitals are from the same I = 2 sub- 
shell, eg, 3d subshell. Since all radial functions are the 
same they factor out from the determinant and for the 
purpose of finding the node they can be omitted. The 
S symmetry allows to rotate the system so that, say, 
electron 1 is on the z-axis, and then the corresponding 
column in the Slater matrix becomes {2z\, 0, 0, 0, 0). As- 
suming that z\ 7^ wc can then write the nodal condition 
as 

det{x 2 — y 2 , xz, yz, xy} = 0. (9) 

Using one of the electrons as a probe (ie, looking at the 
node from the perspective of one of electrons) we can 
find the projection of the node to 3D space. By denoting 
the probe electron coordinates simply as (x, y, z) and by 
expanding the determinant we get 

(x 2 — y 2 )vn\ + xzm 2 + yzm 3 + xyrri^ = (10) 

where rrn are the corresponding cofactors. We divide 
out the first cofactor assuming that it is nonzero (not a 
crucial assumption as clarified below). We get 

(x 2 — y 2 ) + axz + byz + cxy = (11) 

where a — ni 2 jni\, b = m$/m\, c — m^jm\. By com- 
pleting the square this can be further rearranged to 

{x - k lV ){x - k 2 y) + z{ax + by) = (12) 

with fci.2 = (— c± Vc 2 + 4)/2. Let us define rotated and 
rescaled coordinates 

u* = — (afe — b)(x — kiy)/(k\ — fa) (13) 



v* = {ak 1 -b)(x-k 2 y)/{k 1 -k 2 ) (14) 



w* = z[(ak! - b)(ak 2 - b)]/(h - k 2 ) 2 (15) 

so we can write the Eq. fl 1 1 1) as 

u*v* +w*u* +w*v* = 0. (16) 

Note that this equation has a form which is identical to 
Eq. I|1U|) with mi = so this representation is correct for 
general mi. After some effort one finds that Eq. (|16fl is 
a cone equation (ie, d z 2 orbital) as can be easily verified 
by using the following identity 

(2u 2 - v 2 - w 2 )/8 = u*v* + w*u* + w*v* (17) 



where u = u* + v* + 2w* , v — (~u* + v* + 2u>*), w = 
(u* - v* + 2w*). The 3D projected node is therefore 
rotated and rescaled (elliptic) cone. 

At this point it is useful to clarify how the derived node 
projection cone is related to the complete 14-dimensional 
node. Remarkably, the 3D projection enables us to un- 
derstand some of the properties of the 14-dimensional 
manifold. First, the cone orientation and elliptic radii 
(ie, rescaling of the two axes with respect to the third 
one) are determined by the position of the four electrons 
in 3D space: with the exception of special lower dimen- 
sional cases explained below there always exists a unique 
cone given by the Eq. (|16f) which "fits" the positions 
of the four electrons. Besides the special cases (below) 
we can therefore define a projection of a single point in 
4x3= 12-dimcnsional space of four electrons onto a cone. 
That also implies that the complete 12-dimensional space 
describes a set (or family) of cones which are 3D projec- 
tions of the nodal manifold. Similar projection strategies 
are often used in algebraic geometry to classify or an- 
alyze surfaces with complicated topologies and/or high 
dimensionalities . 

Since the cone orientation and two radii are uniquely 
defined by the point in 12 dimensions and the cone it- 
self is a 2D surface in 3D space of the probe electron the 
complete node then has 12+2=14 dimensions. Therefore 
the d 5 HF node is a set of cone surfaces specified by the 
positions of the electrons . This particular form is simply 
a property of the d 5 Hartree-Fock determinant. From the 
derivation above it is clear that after factoring out the ra- 
dial parts one obtains a homogeneous second-order poly- 
nomial in three variables with coefficients determined by 
the positions of the four electrons. In fact, from the the- 
ory of quadratic surfaces^, one finds that a general ellip- 
tic cone can possibly fit up to five 3D points/electrons, 
however, in our case the cone has an additional con- 
straint. Our system was reoriented so that one of the 
electrons lies on the z-axis; that implies that the z-axis 
lies on the cone. Therefore the cone always cuts the xy 
(ie, z — 0) plane in two lines which are orthogonal to 
each other. The orthogonality can be verified by impos- 
ing z — in Eq. (|12f> and checking that k\k 2 = — 1. In 
addition, one can find "degenerate" configurations with 
two pairs of two electrons lying on orthogonal planes (Fig. 
. This corresponds to the " opening" of the cone with 
one of the elliptic radii becoming infinite and the result- 
ing node having a form of two orthogonal planes (Fig. 
|3J). Since in this case there is an additional condition 
on the particle positions, the two-plane node has lower 
dimension and is a zero measure subnode of the gen- 
eral 14-dimensional node. The condition is equivalent 
to A44 = b 2 — a 2 — abc — 0, where A44 is one of the 
quadratic invariants 1 ™. There are more special cases of 
lower dimensional nodes: a) when two electrons lie on 
a straight line going through the origin; b) when three 
electrons lie on a plane going through the origin; c) when 
four electrons lie in a single plane. 

Remarkably, the analysis above enables us to find the 
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number of nodal cells. From Fig. |3|one can infer that by 
appropriate repositioning of the four electrons the cone 
surface smoothly " unwraps" the domains inside the cone, 
forms two crossing planes and then "wraps" around the 
cone domains of the opposite sign. That implies that 
an electron inside one of the cone regions can get to the 
region outside of the cone (with the same wavefunction 
sign) without any node crossing, using only appropriate 
concerted repositioning of the remaining four electrons. 
That enables us to understand that a point in the 15- 
dimensional space (positions of five electrons) can con- 
tinuously scan the plus (or minus) domain of the wave- 
function: there are only two maximal nodal cells. 




FIG. 4: (Color online) Projected Hartree-Fock node of 8 S(/ 7 ) 
state. The node has two topologies: cone times planar surface 
or a cone "fused" with planar surface what forms a single sheet 
surface. There is a smooth transition between these two forms 
depending on the positions of six electrons which are denoted 
by the small spheres. Note that the node contains the z— axis 
which is denoted by the dashed line. 



B. Approximate Hartree-Fock nodes of the 8 S(/ 7 ) 

ion 



planar surface is deformed from a straight plane since it 
passes through the origin and, in addition, it fits three 
of the electrons. The second topology is a "fused" cone 
and planar surface which results in a general single sheet 
cubic surface. The node transforms smoothly between 
these two topologies depending on how the six electrons 
move in space. These two topologies define the projec- 
tion of the node into the probe 3D space and therefore 
enable us to capture the many-dimensional node for this 
particular Hartree-Fock state. This again enables to de- 
scribe the complete node using a theorem from algebraic 
geometry which states that any cubic surface is deter- 
mined by an appropriate mapping of six points in a pro- 
jective planeiiiLiS. To use it we first need to realize the 
following property of the 3D projected node: The node 
equation above contains only a homogeneous polynomial 
in x,y,z which implies that in spherical coordinates the 
radius can be eliminated and the node is dependent only 
on angular variables. Hence, any line defined by an ar- 
bitrary point on the node and the origin (ie, a ray) lies 
on the node. In other words, we see that the surface 
is ruled, ie, it can be created by continuous sweep(s) of 
ray(s) passing through origin. This enables us to project 
the positions of the six electrons on an arbitrary plane 
which does not contain the origin and the node will cut 
such a plane in a cubic curve. As we mentioned above, 
a theorem from the algebraic geometry of cubic surfaces 
and curves says that any cubic surface is fully described 
by six points in a projective plane (seeiii&iS) . For ruled 
surfaces any plane not passing through the origin is a pro- 
jective plane and therefore we can specify a one to one 
correspondence between the 6x3= 18 dimensional space 
and our cubic surface in 3D. Obviously, there will be a 
number of lower-dimensional nodes which will correspond 
to positions of electrons with additional constraints such 
as when they lie on curve with the degree lower than 
cubic; ie, a conic. 



We will use similar strategy as in the preceding case. 
After rotating one of the electrons to z-axis we expand 
the determinant in the probe electron column and elim- 
inate the radial orbitals which form an overall prefactor 
of the Slater determinant since we assume that all seven 
/-states are from the same / = 3 subshell (eg, 4/). We 
get 

(m\x + m2?/)(4z 2 — x 2 — y 2 ) + m^z(x 2 — y 2 ) + m^xyz 



+m 5 x(x 2 - 3y 2 ) + m 6 y{y 2 - 3x 2 ) = (18) 

Note that the node contains the z-axis and there are two 
possible values of z for any x, y since the form is quadratic 
in z. This restricts the node shapes significantly and by 
further analysis one can find that the nodal surface pro- 
jection into 3D has two topologies (Fig. |3}. The first one 
is a cone times a planar surface (topologically equivalent 
to the I40 spherical harmonic). Note that, in general, the 




FIG. 5: (Color online) The 3D projection of the nitrogen 
cation 5 S'(sp 3 ) Hartree-Fock node (the core electrons are elim- 
inated by pseudopotentials) . The projected node exhibits two 
topologies. It is either a planar surface deformed by the radial 
orbital functions at the nucleus or, in certain configurations, 
the deformation forms a small bubble detached form the sur- 
face (the picture on the right). The small cross is the location 
of the ion while the small spheres denote positions of elec- 
trons. 
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C. Approximate Hartree-Fock nodes of the 5 S(sp 3 ) 

ion 



The HF node for this two-shell spin-polarized state can 
be investigated in a similar way as in previous cases with 
a new feature that the radial parts will be present in the 
expansion of the determinant. By expanding the deter- 
minant in the column of the probe electron with position 
x, y, z the 3D node projection is simply given by 



x + b'y + c'z + d'r)(r) = (19) 



where b' 7 c' 7 d' depend on ratios of cofactors and r](r) — 
Ps i r ) I Pp ( r ) is the ratio of radial parts of s and p or- 
bitals and r = \/ x 2 + y 2 + z 2 . The probe electron will 
see a plane with a approximately bell-shape deformation 
in the area of the nucleus (See Fig. [SJ). The shape of 
deformation depends on the ratio of s and p radial parts 
and the magnitudes and signs of the cofactors. For cer- 
tain configurations the deformation is so large that it 
gets detached from the surface and forms a separated 
ellipsoid-like bubble. The bubble results from the ra- 
dial dependence of rj(r) which for pseudized core is not a 
monotonic function and therefore can create new topolo- 
gies. Note that despite the fact that the 3D projection 
shows a separated region of space (the bubble) the com- 
plete node has again the minimal number of nodal cells 
property. To understand this, suppose that the probe 
electron is located inside the bubble and wavefunction 
there has a positive sign. Let us try to imagine how the 
electron can get to the other positive region (the other 
side of the planar surface) . Seemingly, the electron would 
need to cross the nodal surface twice (the surface of the 
bubble and the planar surface). However, the complete 
node is a collective-coordinate object and by moving the 
other two electrons in an appropriate way the bubble at- 
taches to the surface and then fuses into a single surface 
(Fig. 5, left) so that the probe electron can reach the 
positive region without node crossing. 

In order to see whether the correlation would change 
the HF node we have carried out a limited study of the CI 
wavefunction nodes for this case; we have found some dif- 
ferences but we have not discovered any dramatic changes 
to the HF nodes. To quantify this further we have calcu- 
lated the CI energy (with the same basis and level of cor- 
relation as in the previous cases) and the result is in the 
last row of Tab. [I] We estimate that the fixed-node bias 
of the HF node is of the order of ~ 0.001 Hartree which is 
close to « 3 % of the correlation energy. Obviously, the 
DMC energy is above the exact one and percentage- wise 
the amount of missing correlation energy is not insignif- 
icant. We conjecture that the HF node is reasonably 
close to the exact one although the fine details of the 
nodal surface are not captured perfectly. 



D. Approximate Hartree-Fock nodes of 
spin-polarized p 3 d b and sp 3 d 5 shells with S symmetry 

Let us for a moment assume a model wavefunction 
in which the radial parts of s,p,d orbitals are identical. 
Then, using the arrangements similar to d 5 case, we can 
expand the determinant of p 3 d 5 in one column and for 
the 3D node projection we then get 

2u 2 - v 2 - w 2 + au + f3v + -yw = (20) 

where u,v,w are appropriate linear combinations of 
x, y, z. This can be further rewritten as 

2{u + a/4) 2 - (v - (3/2) 2 - (w - 7 /2) 2 + S = (21) 

where S Q = {-a 2 /2 + [3 2 + 7 2 )/4. It is clear that the 
quadratic surface is offset from the origin (nucleus) by 
a vector normal to au + (3v + "/w = plane. Using 
the properties of quadratic surfaces one finds that for 
(a 2 /{a 2 + (3 2 + 7 2 )) < 2/3 the node is a single-sheet hy- 
perboloid with the radius \fSo\ otherwise it has a shape 
of a double-sheet hyperboloid. The double-sheet hyper- 
boloid forms when there is an electron located close to 
the origin. A special case is a cone which corresponds to 
(Sq = 0). The case of sp 3 d 5 is similar, but with different 
5q which now has a contribution from the s-orbital (see 
Fig. EJ). Once we include also the correct radial parts 
of orbitals in the s,p,d channels the coefficients of the 
quadratic form depend on both cofactors and orbital ra- 
dial functions. The resulting nodal surface is deformed 
beyond an ideal quadric and shows some more compli- 
cated structure around the nucleus (see Fig. 0) as il- 
lustrated on HF nodes of the majority spin electrons in 
Mn ++ ion (note that the Ne-core electrons were elimi- 
nated by pseudopotentials) . 



FIG. 6: (Color online) The 3D projection of the angular part 
of the 10 S(sp 3 d 5 ) state Hartree-Fock node (with radial parts 
of orbitals identical for all spd orbitals) . The projection has 
a topology of a single-sheet or double-sheet hyperboloid. The 
small cross shows the location of the nucleus while the spheres 
illustrate the electron positions. 
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FIG. 7: (Color online) Projected Hartree-Fock node of 
10 S(sp 3 d 5 ) of the majority spin valence electrons in Mn +2 
ion. The Ne-core electrons are eliminated by pseudopoten- 
tials. Note the deformations from the radial parts of orbitals, 
including a small bubble detached from the rest of the surface 
(the right picture). For clarity, the positions of electrons have 
been omitted. 



IV. GENERALIZED ELECTRON 
COORDINATES 

What we have learned from the previous cases is that 
for small number of electrons the Hartree-Fock wavefunc- 
tions display nodes which, if transformed in an appropri- 
ate way, lead to rather simple geometries. 

In addition, it is instructive to consider how the nodal 
surfaces evolve with increasing number of electrons. Ob- 
viously, HF theory leads to low kinetic energies and the 
resulting mean-field nodes are very smooth. The exact 
nodes of the high symmetry P{p 2 ) and S(p 3 ) cases can 
be interpreted as reoriented planes which enable us to 
fit one or two electrons and obviously such rotations do 
not cause any increase in the kinetic energy. For more 
particles the rotations and translations are not sufficient 
to fit the electron positions and the lowest increase in 
kinetic energy in atomic systems is apparently produced 
by rescaling of the axes. 

Finally, for larger number of electrons the node be- 
comes more deformed and complex with possibilities of 
new topologies and topological changes. The 3D node 
projections we have analyzed above show that often 
there exist coordinate transformations which can sim- 
plify the node description and enable us to find useful 
node parametrizations (at least, for our cases of spin- 
polarized electronic subshells in ions). As we demon- 
strated on the 3 P(p 2 ) state similar spin-polarized open 
shells can be studied by analogous techniques as well. 

In the analysis of the nodes we have presented a num- 
ber of coordinate transformations and maps which en- 
abled us to understand particular nodal structures and 
their topologies. It is interesting to explore this idea fur- 
ther and think about possible research directions. In or- 
der to illustrate some of the possibilities let us define a 
single-electron coordinate — > r* transformation as 

r* = M i (R)r i + t i (R) (22) 



where Mj(i?) is a metric tensor, tj(i?) is a spatial offset 
(translation) while R represents coordinates of all elec- 
trons. The dependence of Mj(iZ) and tj(i?) on R can 
be nonlinear and even include an explicit dependence on 
external potential to describe charge inhomogeneities or 
required symmetries. The antisymmetry condition furher 
restricts the dependences of Mj(i2) and t,(ii) on R. Ob- 
viously the metric tensor has to be positive definite but 
otherwise the transformation is variationally free. Actu- 
ally, the usefulness of this transformation can be twofold. 
First, as we mentioned, it can be employed as an efficient 
way to project the high-dimensional nodal manifold into 
a simpler, low-dimensional projection. Second, it can 
be useful as a new variational freedom to optimize the 
wavefunction nodes using the following form of a single 
determinant or linear combination of determinants 

* = J2dkdet k {^(r*)} (23) 

k 

where {<Pa^} are one-particle orbitals and dk are expan- 
sion coefficients. Therefore, besides optimizing the or- 
bitals, one can also optimize the metric tensor and off- 
set in order to get better variational wavefunctions and 
fixed- node energies. In fact, the transformation above 
can be considered a generalization of the Feynman-Cohen 
backflow quasi-particle coordinates^. By simplifying 
Mj(.R) to the unit matrix times a scalar function we 
can easily recover the backflow wavefunction which Feyn- 
man and Cohen suggested for liquid helium 23 and which 
was successfully employed in QMC in several previous 
studies 2 ^ 2 ^^. The new feature proposed here is the met- 
ric tensor which enables to better describe the systems 
with inhomogeneities and/or with rotation symmetries. 



V. CONCLUSIONS 

We have investigated the nodes of atomic and molec- 
ular spin-polarized systems with one-particle states in 
s,p,d channels. We have studied cases with high sym- 
metries which enabled us to find exact nodes for sev- 
eral states with a few electrons (p 2 ,p 3 ,ir 2 ). Moreover, 
the projection of multi-dimensional manifolds into 3D 
space enabled us to study and characterize properties 
of nodes, in particular, their topologies for the Hartree- 
Fock wavefunctions. This analysis has provided useful in- 
sights and enabled us to formulate a general transforma- 
tion of one-particle coordinates using coordinate trans- 
lation (backflow) and metric tensor to capture inhomo- 
geneities and/or rotation symmetries. Such transforma- 
tions can be useful for understanding the nodal properties 
and topologies and also as a new variational freedom for 
QMC trial wavefunctions. 
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